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ABSTRACT 

Viscous redistribution of mass in Milky Way-type galactic disks is an appealing way of generating an expo- 
nential stellar profile over many scale lengths, almost independent of initial conditions, requiring only that the 
viscous timescale and star formation timescale are approximately equal. However, galaxies with solid-body ro- 
tation curves cannot undergo viscous evolution. Low surface brightness (LSB) galaxies have exponential surface 
brightness profiles, yet have slowly rising, nearly solid-body, rotation curves. Because of this, viscous evolution 
may be inefficient in LSB galaxies: the exponential profiles, instead would give important insight into initial 
conditions for galaxy disk formation. 

Using star formation laws from the literature, and tuning the efficiency of viscous processes to reproduce an 
exponential stellar profile in Milky Way-type galaxies, I test the role of viscous evolution in LSB galaxies. Under 
the conservative and not unreasonable condition that LSB galaxies are gravitationally unstable for at least a part 
of their lives, I find that it is impossible to rule out a significant role for viscous evolution. This type of model still 
offers an attractive way of producing exponential disks, even in LSB galaxies with slowly-rising rotation curves. 

Subject headings: galaxies: general — galaxies: evolution — galaxies: structure — galaxies: spiral 



1. INTRODUCTION 

An exponential stellar light profile over 4-6 disk scale lengths 
is an almost universal observational feature of disk galaxies 
' (e.g.. Freeman 1970; de Jong 1996). Yet, the production of 
. a stellar exponential disk over more than 3 scale-lengths dur- 
' ing the early stages of galaxy formation can prove highly chal- 
lenging (e.g., Dalcanton, Spergel & Summers 1997; van den 
Bosch 2001; Ferguson & Clarke 2001, although see also Con- 
tardo, Steinmetz & Fritze-von Alvensleben 1998). Ever since 
their inception (Lin & Pringle 1987a), viscous evolution mod- 
els have endured as an attractive way of producing exponential 
stellar disks over >4 disk scale lengths, almost independent of 
. the initial density distribution. In a differentially rotating gas 
' disk, viscosity caused by non-circular gas motions and turbu- 
lence transports angular momentum outwards as mass flows in- 
wards (Pringle 1981; Lin & Pringle 1987a). To prevent viscous 
evolution from reaching its logical endpoint (all the mass at the 
origin, all the angular momentum at infinity), the star forma- 
tion (SF) timescale t^, should be within half an order of magni- 
tude of the viscous timescale' ~ /v to 'freeze in' a nearly 
exponential stellar surface brightness profile over many disk 
scale-lengths (e.g., Lin & Pringle 1987a; Clarke 1989; Yoshii & 
Sommer-Larsen 1989; Olivier et al. 1991; Hellsten & Sommer- 
Larsen 1992; Firmani, Hernandez & Gallagher 1996; Ferguson 
& Clarke 2001; Slyz et al. 2002). 

In ID, the evolution of the gas density Sg at a given radius r 
is given by: 



dUg _ Id Udldr)[vY.gr^dn/dr] 
~W~~r^r\ d/dr{r^n) 




where t is the time, v is the viscosity, is the angular velocity, 
and "0* is the SF rate (Lin & Pringle 1987a). A brief inspection 
of Eqn. 1 shows a well-known result: in a galaxy without shear 
(a solid-body rotation curve with dil/dr = 0) there can be no 
viscous evolution. 

' This expression for only applies in the case of a rotation curve which is 
not close to solid body. 



Bearing this in mind, it is interesting to note that recently 
published high-resolution (<2") Ha rotation curves for low 
surface brightness (LSB) galaxies are slowly rising, and in a 
few cases are nearly solid-body, within the optical extent of 
the galaxy (e.g., de Blok et al. 2001; de Blok & Bosma 2002; 
Matthews & Gallagher 2002). Solid-body rotation curves have 
no shear, meaning that viscous evolution would not occur Yet, 
many LSB galaxies have exponential surface brightness profiles 
over 4 or more disk scale lengths (e.g., McGaugh & Bothun 
1994; de Blok, van der Hulst & Bothun 1995; Bell et al. 2000). 
This could argue against a significant role for viscous evolution 
in generating the exponential stellar disks of LSB galaxies, im- 
plying that the structure of (at least) LSB disks offers important 
insight into the initial conditions of galactic disk formation. 

To understand if constraints on galaxy formation models can 
be gleaned from the structure of LSB galaxy disks, it is impor- 
tant to quantitatively explore this issue. The key questions are: 
i) how close to solid-body are the rotation curves of LSB galax- 
ies? and ii) will the same kind of viscous evolution processes 
postulated to affect Milky Way-type galaxies be able to capital- 
ize on any small departures from solid-body behavior, and be 
able to generate exponential stellar disks in LSB galaxies? 

To address this problem, I set out in this paper to estimate 
the possible effect of viscous evolution in LSB galaxies with 
known, well-constrained surface brightness profiles and rota- 
tion curves. I make the central assumption that viscous pro- 
cesses in LSB galaxies and normal Milky Way-like galaxies 
have the same efficiency, and that any differences in viscous 
timescales are the result of different rotation curves and gas 
densities only. Comparing these viscous timescales for LSB 
galaxies with realistic SF timescales (e.g.. Bell et al. 2000), one 
can say whether the same sources of viscosity that are applica- 
ble in normal galaxies could affect the evolution of LSB galax- 
ies. 

The plan of this paper is as follows. In §2, a brief description 
of the model is given. In §3, the model is caUbrated to pro- 
duce an exponential stellar light distribution for the Milky Way. 
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In §4, the Milky Way-calibrated model is applied to five LSB 
galaxies. I present my conclusions in §5. As an aside, I discuss 
viscosity from gas cloud collisions in Appendix A. 

2. THE MODEL 

To investigate the role of viscous evolution in LSB galaxies, 
I modify the model of Bell & Bower (2000) to include viscous 
evolution. Following Lin & Pringle (1987a), I adopt a smooth 
initial gas surface density distribution Sg(r): 

^ r So[l+cos(7rr/5/i)] 0<r<5h 

1^ r> 5h, ' 

where r is the radius, h approximates the final exponential 
scale-length, and 2S() is the central surface density (the final ex- 
ponential central surface density would be roughly SEq; see the 
dashed line in Fig. 1). The choice of this particular initial pro- 
file is arbitrary: Yoshii & Sommer-Larsen (1989) demonstrate 
that exponential disks are produced by viscous evolution from 
any smooth, centrally-concentrated initial profile. The evolu- 
tion of gas density with time as a function of radius is followed 
using Eqn. 1. 

SF is followed using either i) a gas density-dependent SF 
law (e.g., Schmidt 1959) V* = kY."^, where k is the rate of SF 
at a gas surface density of 1 Mqpc"^, and n dictates how 
sensitively SF rate depends on gas surface density, or a 
dynamical timescale-dependent SF law (e.g., Kennicutt 1998) 
^* = KT,g/Tdyn, where K is the rate of SF at a gas surface den- 
sity of 1 M0 pc~^ and dynamical time of 1 Gyr, and Tdyn = 
6.16r(kpc)/y(kms~ ) Gyr is the time taken to orbit the galaxy 
at a distance r. Both SF laws are reasonably consistent with the 
SF rates and histories of present-day spiral galaxies (e.g., Ken- 
nicutt 1998; Bell & Bower 2000). In the following, values of 
k, n, K and the chemical element yield are fixed at the values 
quoted in Table 1 of Bell & Bower (2000). 

At present, sources of viscosity in galactic disks are not well- 
understood. In previous studies, the SF and viscous timescales, 
f* and t,j, were directly related (where, again, only applies 
where the rotation curve is not nearly solid body). The propor- 
tionality constant (i = t^,/ti, was constrained to be within half 
an order of magnitude of unity in order to produce exponen- 
tial stellar disks over 4-6 disk scale lengths (e.g., Lin & Pringle 
1987a; Hellsten & Sommer-Larsen 1992; Ferguson & Clarke 
2001). In this work, I wish to relax this constraint to a certain 
extent. In particular, one of the things that I wish to assess is if 
(3 for LSB galaxies, as this would imply a minor role for 
viscous evolution in LSB galaxies. 

In this paper, I explore viscosity due to gravitational instabil- 
ities. Since I find later that viscous evolution in Milky Way-type 
and LSB galaxies is efficient with this physically-motivated 
prescription, it is sufficient to adopt only one viscosity prescrip- 
tion for the purposes of this paper. Following, e.g., Olivier et 
al. (1991) 1 also explored viscosity due to collisions between 
gas clouds; however, I found that cloud-cloud colhsions are in- 
capable of driving significant viscous evolution in any spiral 
galaxies (see Appendix A for the derivation of this result). I 
therefore do not consider this prescription further. Large-scale 
gravitational instabilities give a viscosity (Lin & Pringle 1987b; 
Olivier et al. 1991): v = agrav7r^G^E^/K^ where agrav is a con- 
stant, G is the gravitational constant, and Eg is the gas density 
in M© pc~^. The epicychc frequency k in Gyr~^ is given by 
K = {r d / dr{^^)+A^^YI''- . I set the viscosity from gravitational 
instabilities to zero if the disk is gravitationaUy stable, i.e., if the 



Toomre (1964) Q> I, where Q = VsK/nCSg. The gas velocity 
dispersion Vj is assumed to be 6kms"^ (Kennicutt 1989; Bin- 
ney & Tremaine 1987). The constant agrav is set to produce an 
exponential stellar disk for the Milky Way (essentially, /? for 
the cloud-collision and gravitational instability cases is set to 
be ~ 1: see §3). This value of agrav is then left alone, and the 
LSB galaxy models are run to test if the difference in rotation 
curve shape and gas density significantly affects (3. This is the 
central assumption of this analysis. 

Stellar populations are modeled using the Bruzual & Chariot 
(in preparation) stellar population synthesis models (see, e.g., 
Liu, Chariot & Graham 2000) adopting a Salpeter IMF with a 
reduced number of low mass stars (following Bell & de Jong 
2001). 1 adopt the instantaneous recycling approximation and 
the closed box approximation (Bell & Bower 2000). The de- 
tailed viscous flow of metals is not tracked in this model. These 
approximations do not affect the A'-band surface brightness or 
B-R color significantly, and therefore will not affect the con- 
clusions of this paper. It is worth noting that use of the B-band 
or stellar mass surface density profiles would also not affect the 
conclusions of this paper 

1 solve these equations using a standard first order explicit 
scheme with 50 equally-spaced radius steps between 0.2h and 
IQh. Time steps are 4 Myr, for a galaxy age of 12 Gyr. There are 
two runs per galaxy model with the different SF laws. Follow- 
ing e.g., Lin & Pringle (1987a) and Ferguson & Clarke (2001), 
1 allow mass to be lost from the center of the galaxy (i.e., total 
galaxy mass is not conserved); however, in practice less than 
10% of the total mass is lost for even large amounts of viscous 
evolution. This does not significantly affect my conclusions, as 
the galaxies with large amounts of viscous evolution are typi- 
cally quite luminous, and this 'lost' mass is easily accommo- 
dated within a bulge component. 

3 . REPRODUCING THE MILKY WAY 

For the Milky Way, I adopt a flat rotation curve at all radii 
with V = 220kms~^ (adopting a more realistic rotation curve 
in the inner parts does not significantly change the results). I 
choose a value of h = 2.82 kpc, and Eq = 200 Mq pc"^, resulting 
in disk central surface densities of ~ lOOOM© pc~^ (see Table 
1 for the modeUng parameters of this and subsequent galaxy 
models). These values of V, h and the disk central surface den- 
sity are in reasonable accord with observations (e.g., Binney 
& Merrifield 1998). 1 assume that the gas disk was assembled 
instantaneously 12 Gyr ago (assuming more realistic infall his- 
tories for any of the models in this paper does not significantly 
change the results). A choice of Agrav = 0.02 for the gravita- 
tional instability model results in a profile which is exponential 
to within 0.11 mag at r < 5/z for the gas-density dependent SF 
law, and to within 0.14 mag for the dynamical time SF law. 
In bofli cases, B-R = 1.1 ±0.02, tiie gas fraction is 23% ± 
3%, and the /iT-band central surface brightness is 17.35 ±0.15 
mag arcsec"^. These global parameters are typical of a galaxy 
of the Milky Way's luminosity and size (e.g.. Bell & de Jong 
2000). Variations of Ggrav by factors of 2-3 do not significantly 
affect the surface brightness profiles of the model Milky Way- 
type galaxies: they are still exponential to <20% within 4 disk 
scale lengths. 

There are two interesting points which deserve further dis- 
cussion. Firstly, my value of agrav results in much slower vis- 
cous evolution than estimated by Olivier et al. (1991). Ohvier 
et al. (1991) estimate short viscous timescales <0.2 Gyr for gas 
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Table 1 

Parameters for galaxy models 





En 


h 


V 

(kms"') 






Milky Way 


200 


2.82 


220 






ESO-LV 1870510 


16 


3.46 


52(l-e-'-[kpc]/2.8) 






UGC 11557 


120 


2.77 


110(l-e-'^P'^l/5') 






F568-3 (soUd-body) 


40 


2.96 


19.95 r[kpc] 




r< 3.85kpc 








76.8 + 9.22(r[kpc]- 


-3.85)1/2 


r > 3.85kpc 


F568-3 {la curvature) 


40 


2.96 


1 9.95 rfkpcl + 6(1 - 


0.27{r[kpc]- 1.923 }2) 


r < 3.85kpc 








76.8+9.22(r[kpc]- 


-3.85)'/2 


r > 3.85kpc 


F583-1 


20 


1.71 


90(l-e-'-[kpc]/3-26) 






F583-4 


16 


2.14 
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Fig. 1. — AT-band surface brightness profile of a Milky Way model with 
a gas density-dependent SF law, which evolves subject to viscous evolution 
from gravitational instabilities. The solid fine denotes the final S'-band sur- 
face brightness profile, the dotted line denotes an exponential fit to the sur- 
face brightness profile for r < 6/1, the dot-dashed line shows the star-Fgas sur- 
face density profile assuming a /("-band stellar M/L of 1 in solar units, and the 
dashed line denotes the initial gas density profile, again assuming M/Ljs: = 1 . 
The dynamical time SF law yields similar results. 



at '^Skpc from the galactic center. In order to keep l3 ^ I, 
it is necessary to have SF timescales ^ 0.2 Gyr. These short 
timescales violate measured SF timescales in large spiral galax- 
ies in general (e.g., Kennicutt 1998; Bell & de Jong 2000), and 
in the Milky Way in particular (> 10 Gyr in the solar cylin- 
der; Rocha-Pinto et al. 2000). In addition, their viscous flows 
would have velocities ^ 5 km s"' , violating the < 1 km s~' con- 
straint derived from radial metallicity distributions (Lacey & 
Fall 1985; Clarke 1989). In stark contrast, the 'low' viscous 
efficiencies that I adopt result in viscous and SF timescales of 
~ 10 Gyr at the solar cylinder, and gaseous flow velocities of 
~0.1kms"i: well within observational constraints. Secondly, 
the model surface brightness profile is exponential to ~ 1 0% de- 
spite the ratio of the SF and viscous timescales, /? = f*/f,y, vary- 
ing with radius by nearly an order of magnitude. This strength- 
ens the conclusion of HeUsten & Sonimer-Larsen (1992), who 



stated that variations in /3 by a factor of two with radius resulted 
in an exponential profile: in fact, an order of magnitude varia- 
tion in (3 with radius is acceptable in some situations. 

4. ARE LSB GALAXIES AFFECTED BY VISCOUS EVOLUTION? 

The stage is now set to address the role of viscous evolution 
in LSB galaxies. All of the ingredients are in place: i) accurate 
observations of rotation curves for real LSB galaxies with well- 
constrained exponential profiles (required for Equation 1), /;) 
knowledge about SF timescales in LSB galaxies and plausible 
SF laws (Bell & Bower 2000), and in) estimates of the effi- 
ciency of gravitational instability viscosity in Milky Way-type 
galaxies from §3. The central assumption of this comparison is 
that the efficiency of viscous processes in LSB galaxies is the 
same as for Milky Way-type galaxies (i.e., agrav is the same for 
all galaxies). With this assumption, I am postulating that the 
physical processes in LSB and Milky Way-type galaxies are 
identical, and differences in behavior result from differences in 
rotation curve and gas density only. 1 will present the anal- 
ysis for one particular LSB galaxy, ESO-LV 1870510, before 
extending the analysis to a further four suitable LSB galaxies. 

4.1. ESO-LV 1870510 

ESO-LV 1870510, at a distance of 29 Mpc (Bell et al. 2000), 
has a slowly-rising, almost solid-body, Hci-derived rotation 
curve from de Blok et al. (2001, see their Fig. 1 and the in- 
set panel in Fig. 2). Bell et al. (2000) present optical and 
near-IR photometry (their Fig. 1), which demonstrates that 
the galaxy is exponential to within 0.1 mag over 4 disk scale 
lengths. ESO-LV 1870510 has an intrinsic B-band central sur- 
face brightness of 23 .75 ±0.2 mag arcsec"'^, implying nearly a 
factor of 10 lower surface density than the canonical 21.65 
magarcsec^- of Freeman (1970). For the model, I adopt the 
following smooth approximation to de Blok's smoothed rota- 
tion curve: V = 52(1 - e~''^^'^^y^'^)kms~^ (see inset to Fig. 2), a 
value of So = 16 M© pc"'^ and h = 3.46 kpc (see also Table 1). 

Fig. 2 shows the resulting A'-band surface brightness pro- 
files of the density-dependent and dynamical time SF law LSB 
model galaxies. Clearly, the viscous redistribution and SF 
timescales turn out to be roughly equal, and exponential surface 
brightness profiles are produced. The stellar disks are exponen- 
tial to within 0.07 mag (density-dependent SF law) or 0.13 mag 
(dynamical time SF law) within r < 4h. The model B-R col- 
ors match the observations to < 0. 1 mag, and the scale lengths, 
magnitudes and central surface brightnesses match the observa- 
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Fig. 2.— A:-band surface brightness profiles of two ESO-LV 1870510 mod- 
els which evolve subject to gravitational instability viscosity. The gas density- 
dependent SF law is shown as a thick solid line and the dynamical time SF 
law is shown by a thick dashed line. The exponential fits to both models are 
also shown (thin soUd and dashed lines). The dotted line denotes the initial gas 
density profile, assuming M/Lk = 1. Inset is the smoothed rotation curve of de 
Blok et al. (2001) and my fit to the rotation curve (solid line). 



tions to 20%. ESO-LV 1870510 lacks literature Hi data, how- 
ever, the model gas fraction of ~60% is typical of LSB galaxies 
(e.g.. Bell & de Jong 2000). 

4.1.1. Why Can Viscous Evolution Affect LSB Galaxies? 

It is interesting to consider why viscous evolution can work 
with similar efficiency in LSB galaxies and Milky Way-type 
galaxies. After all, the rotation curves of LSB galaxies rise 
slowly (see, e.g., the inset panel in Fig. 2), implying that viscous 
evolution could be much less efficient. Using that ti, ~ r^/i^, 

oc r^K^/S^,. Fig. 3 of de Jong & Lacey (2000) shows that 
Milky Way-type and LSB galaxies have similar scale lengths, 
therefore characteristic values of r are similar for LSB and 
Milky Way-type galaxies (see also Table 1). The gas density 
Y,g is a factor of ^ 10 lower in LSB galaxies, which drives up 
tv However, k ^ oc V/r (Binney & Tremaine 1987) is sub- 
stantially lower in LSB galaxies as r is similar but V is a factor 
of a few lower. Consequently k^/S^, therefore t^, are similar 
for Milky Way-type and LSB galaxies. 

Considering the rest of Equation 1, we see that the top term 
in the equation is reduced for a LSB galaxy by a factor of 100: 
a factor of 10 from the term and a factor of 10 from dVl/dr. 
The bottom term is reduced by a factor of '--^5, simply due to 
the more solid-body rotation curve shape. Therefore the ratio, 
and the derivative of the ratio, is roughly a factor of 20 lower 
for a LSB galaxy. However, the gas densities are roughly a 
factor of 10 lower, so {dTig/ dt} /Tig is only reduced by a factor 
of 2 or so: viscous evolution from gravitational instability will 
be nearly as effective in LSB galaxies as it would be in large 
spiral galaxies! Coupled with the somewhat longer in LSB 
galaxies, it is clear that (3 = ti,/tv will remain of order unity, 
easily producing an exponential disk. 



4. 1 .2. Gravitational Instability and LSB Galaxies 

I have demonstrated that viscous evolution through gravita- 
tional instability can be just as effective in LSB galaxies as 
in Milky Way-type galaxies. However, this seems to contra- 
dict at some level the argument that LSB galaxies are rather 
more likely the be gravitationally stable than their higher sur- 
face brightness counterparts. If LSB galaxies were gravitation- 
ally stable, no viscous evolution (at least through gravitational 
instabihties) would occur, van der Hulst et al. (1993) estimate 
the Toomre Q parameter, finding 1 < 2 ^ 3. This can be com- 
pared to 0.5 < 6 ^ 2 for most luminous spiral galaxies in their 
inner parts (Kennicutt 1989). However, there are reasons to 
believe that van der Hulst et al. (1993) could significantly over- 
estimate Q. Toomre's Qocl /T,g, and poor resolution Hi obser- 
vations will tend to underestimate Eg. In addition, it is unclear 
how much molecular hydrogen there is in LSB galaxies, due 
primarily to the uncertainty in the CO to H2 ratio X (see, e.g., 
Mihos, Spaans & McGaugh 1999; Boselli, Lequeux & Gavazzi 
2002): a 50% increase in Eg caused by molecular hydrogen 
in the inner few scale lengths is not unreasonable given that 
X could be a factor of 10 higher than galactic, and given the 
level of current CO detections of LSB galaxies (Matthews & 
Gao 2001). Both of these effects would tend to drive down Q, 
making LSB galaxies less stable. 

However, the model galaxy has gas densities which are some- 
what larger than typical observed densities in LSB galaxies 
(10 M©pc"^ at the model half-hght radius, compared to ~5 
Mq pc~^ for observed LSB galaxies). This problem can be alle- 
viated by allowing the model galaxy to be built up by a lengthy 
gas infall. Assuming an 'inside-out' formation scenario with an 
e-folding timescale of 2.5 Gyr at the center, and 10 Gyr at ~ 4h 
(the infall case of Bell & Bower 2000), the gas density of the 
model galaxy is reduced by a factor of two, bringing it much 
closer to typical, observed gas densities in LSB galaxies. How- 
ever, viscous evolution still occurs, producing a stellar expo- 
nential disk: the viscous evolution basically operates at a level 
which keeps the LSB galaxy on the edge of gravitational insta- 
bility throughout its life. Given the considerable observational 
and modeling uncertainties, it is not unreasonably conservative 
to argue that LSB galaxies are at the edge of gravitational in- 
stability for most of their lives (note that Mihos, McGaugh & 
de Blok 1997, also argue for locally unstable LSB disks), and 
therefore could plausibly undergo significant viscous evolution. 

4.2. Other LSB galaxies 

I have demonstrated that it is possible for the LSB galaxy 
ESO-LV 1870510 to undergo significant viscous evolution if it 
is gravitationally unstable for at least part of its life. In this 
subsection, I repeat the analysis with a further four galaxies 
with accurate, modelled rotation curves (de Blok et al. 2001) 
and exponential profiles to good accuracy over at least three 
disk scale lengths (de Blok, McGaugh & van der Hulst 1996; 
Swaters 1999; Bell et al. 2000): UGC 11557, F568-3, F583-1 
and F583-4 (see Table 1 for modeling parameters and Fig. 3 for 
the results). 

In all cases the gas disks can be unstable, and this leads to sig- 
nificant viscous evolution with either the gas density-dependent 
or dynamical time-dependent SF law. The cloud collision vis- 
cosity model gives very little viscous evolution. None of the 
rotation curves are solid body: most rise slowly out to ^2h, and 
turn over at larger radii. This is common amongst LSB galaxies 
(e.g., de Blok et al. 2001; Matthews & Gallagher 2002). The 
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Fig. 3. — A^-band surface brightness profiles of the models of four LSB 
galaxies from the sample of de Blok et al. (2001). The gas density-dependent 
models, which evolve due to viscosity from gravitational instabilities, are 
shown for all galaxies. F568-3 has two models shown: one with a solid-body 
inner rotation curve (solid line) and one with curvature which is consistent 
with the Icr errors in the rotation curve (dashed line). With the exception of 
the soUd-body case of F568-3, all surface brightness profiles do not deviate 
from exponential by more than ~0.1 mag. The use of the dynamical time SF 
law would not affect the results significantly. 



only exception to this is F568-3, which has a rotation curve 
within ~1.2/i which can be modelled as solid body. This leads 
to a 'ring' feature in the viscous model surface brightness pro- 
file, where material piles up at r ^ /z (Fig. 3, solid line with 
stars). If the rotation curve is allowed to 'curve' slightly at 
the observational la level, the final surface brightness profile 
smooths considerably, and becomes exponential to within 0.1 
mag RMS (Fig. 3, dashed line with stars). 

This strengthens the conclusions considerably: it is impossi- 
ble to rule out a significant role for viscous evolution during the 
life of five LSB galaxies with exponential surface brightness 
profiles for which accurate rotation curves are available from 
the literature. 

5. CONCLUSIONS 

Models of galaxy evolution which include viscosity due to 
non-circular and turbulent gas motions have, for the last 15 
years, provided a relatively appealing way of producing accu- 
rately exponential stellar disks in large spiral galaxies almost 
independent of initial conditions. The primary assumption is 
that the viscous timescale and star formation timescale should 
be approximately equal. However, galaxies with soUd-body 
rotation curves cannot undergo viscous evolution. Low sur- 
face brightness (LSB) galaxies have exponential surface bright- 
ness profiles, yet have slowly rising, nearly solid-body, rotation 
curves. Thus, it could be that LSB galaxies have had little redis- 
tribution of gas due to viscous evolution during their lifetimes. 
Their surface brightness profiles would give unique insight into 
the initial conditions of galaxy formation. 

In this paper, I have constrained the importance of viscous 
evolution during the life of LSB galaxies. I have adopted star 



formation laws from the literature, and have tuned the efficiency 
of viscous processes to reproduce an exponential stellar profile 
in Milky Way-type galaxies. Conservatively, I have assumed 
that the physics of turbulence are the same in LSB and Milky 
Way-type galaxies: any differences in the importance of vis- 
cous evolution in the two types of galaxy are the product of 
differences in rotation curve and gas density only. I then apply 
this model to LSB galaxies with observed rotation curves and 
surface brightness profiles. Under the conservative and not un- 
reasonable condition that LSB galaxies are gravitationally un- 
stable for at least part of their lives, 1 find that viscous evolu- 
tion from gravitational instabilities is quite effective. Viscous 
evolution models still offer an attractive way of producing an 
exponential disk over four or more disk scale lengths, almost 
independent of initial density distribution, even for LSB galax- 
ies with slowly-rising rotation curves. 

Does this have implications for the power of surface bright- 
ness profiles to constrain galaxy formation theory? It is impos- 
sible at present to predict the degree of viscous evolution that 
will occur, as the physical mechanisms driving viscous evolu- 
tion are so poorly constrained. However, this paper has shown 
that it can plausibly occur in both low and high surface den- 
sity disk galaxies, and Yoshii & Sommer-Larsen (1989) have 
shown that it would modify most intial density profiles into 
an exponential disk. This reduces substantially the ability of 
observations to rule out models of galaxy formation based on 
surface brightness profiles. Indeed, the assembly of a disk 
with a reasonably centrally concentrated profile with roughly 
the right amount of angular momentum is all that is required 
of a galaxy formation model. However, it is impossible to re- 
duce the central matter density of a galaxy through viscous pro- 
cesses; therefore, models which overpredict the central mass 
density of galaxy disks cannot escape the need for revision or 
modification (see, e.g., van den Bosch 2001). 

I thank Dennis Zaritsky, Erwin de Blok and Stephane Court- 
eau for useful discussions during the course of this work. Spe- 
cial thanks go to the anonymous referee, whose insightful com- 
ments led me to examine cloud-collision viscosity more closely. 
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APPENDIX 

A. EXPLORING CLOUD-CLOUD VISCOSITY 

Here, I briefly discuss cloud-cloud collision viscosity. De- 
spite its frequent use in the literature, it has not been com- 
prehensively derived and so I will explore it here. Following, 
e.g.. Silk & Norman (1981) and Olivier et al. (1991) I assume 
a cloud-cloud collision rate R ^ «cioiidKs7rr^iou(j, where Wcioud is 
the number space density of clouds, is the gas velocity dis- 
persion, and Tcioud is the typical cloud radius. Assuming that 
the cloud properties and galaxy scale height are independent of 
gas surface density, ndoud = ^g/(Mc\oudh), where E^, is the gas 
surface density, Mdoud is the typical gas cloud mass and h is the 
typical gas scale height. Adopting the cloud-colUsion parame- 
ters from OUvier et al. (1991), R ~ 0.032vsEg Gyr^i where 
is in Mq pc"^ and v, is in km s"' (so far, this derivation does not 
significantly deviate from Olivier et al. 1991). 

Following Pringle (1981), the kinematic viscosity takes the 
basic form u v,A, where A is the effective mean free path 
of clouds before they collide. In the limit that cloud collisions 
are much more frequent than once per orbit, A is equal to the 
mean free path of the clouds which is VjTc, where = l/R 
is the cloud-collision timescale. This yields a viscosity = 
32.6r7freqi'.s/Sgkpc^Gyr"'. In the limit of cloud collisions be- 
ing much less frequent than once per orbit, v ^ v^XriTc/TK)'^ 
where Ar ~ v?Tk is the radial excursion during one epicycle, 
and Tk = 27r/K is the epicycle timescale (Pringle 1981). In this 
case, V ~ v^t^/tc = 1.3l77infreqVjSg/K^kpc^Gyr"\ where Vj is 
in kms~', is in Mq pc"^, and k is in Gyr"'. Which limit is 
adopted depends on Eg (assuming that is constant between 
different galaxy types). For gas densities >25 Mq pc"^, the fre- 
quent collision case is more appropriate (t^ < 200 Myr), and for 
gas densities <25 Mq pc~^ the infrequent case is more appro- 
priate (tc > 200 Myr), for a constant i\ ^ 6 km s"' . 

The crucial issue now is to estimate the efficiencies r^fteq 
and ryinfreq of the frequent and infrequent cloud-cloud collision 
cases. This is done by estimating the energy loss per cloud- 
cloud collision. The specific kinetic energy of a cloud is ^ v^q,, 
where Vrot is the rotation velocity. However, clouds collide with 
speeds much less than than Vrot (because they are almost al- 
ways in very similar, nearly circular orbits). In this case, the 
maximum possible specific energy loss is related to their rela- 
tive speeds ~ v^^i^^^e- Therefore, the efficiency of energy loss 
per collision is (Vreiative/vrot)^- The relative velocity Vreiative 



can be approximated by the difference in velocity caused by 
the radial excursion of the cloud. Using that Vreiative = VrotAr/r, 
where r is the radius, and substituting v^t = ri7 (O is the an- 
gular frequency) and Ar ~ A, Vreiative ~ f^A. For the frequent 
collision case A ^ TcV,,, therefore Vreiative ^ ri/0.032Eg (where 
17 is in Gyr-i and E^ is in Mq pc'^). Thus, rjfreq ~ 911 /{r^'T?-), 
where r is in kpc. For the infrequent collisions case A is the 
maximum epicyclic excursion, which is roughly v, multiplied 
by the time for 1/4 epicycle ~ Vj/k. Using that ^ 1, 
^relative ^ for the infrequent case, where v, is in km s~^ . There- 
fore, 7/i„fteq ~ V^/(r^r2^). 

What are reasonable estimates of the timescales of viscous 
evolution in the frequent (Milky Way-type galaxy) and infre- 
quent (LSB galaxy) cases? Using r = 1.5 kpc, Vrot = 220 km s"\ 
V, = 6 km s~^ , and = 50 Mq pc"^ for a Milky Way-type galaxy, 
I find that v,-eiative ^ 20kms"' and f,^ ^ r^/i^ ^ 2000 Gyr Sim- 
ilarly, adopting E^ = IOMqPC"^ and Vrot = lOOkms"' for the 
LSB galaxy infrequent case, I find ti, ~ 1000 Gyr. Thus, it is fair 
to conclude that viscosity from cloud-cloud collisions will be 
unimportant for all spiral galaxies. This is because the low effi- 
ciency of energy loss in cloud-cloud collisions (vreiative/ t^'rot)^ ^ 
0.01, which more than outweighs the fact that cloud collisions 
tend to happen > 1 time per orbit for most spiral galaxies. 
Of course, if some mechanism 'whips up' cloud collision ve- 
locities (such as a galaxy interaction), it may be possible to 
evolve significantly from cloud-cloud viscosity, but this mecha- 
nism should not play a significant role in the evolution of spiral 
galaxies in their normal, quiescent mode. 



